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ABSTRACT 

‘r figur “ i0 " s 10 “ ntroI latoge. 

recognition of a deficit? the exisfin^nm^ri 001110 * , u*a r 0 ' 0 ^"™^ response. In 
effort was established in 1990 by NAsK A° l0gy f ° r SeaIs, a seven y ear 

Technology, under the Earth tn nrhi? p m « of Aeron audcs Exploration and 
Computational Fluid Dynamics ^ program to develop validated 

provide NASA and the U S AeSL?Ssi^ w?th^ anaIy i es ^L seaIs - ' n * effort will 
industrial codes for analyzing , CFD “ K " Uflc “*• 

^^ V tat!oMl 3 meSoIogy ^ forflowln^si^ t vdoped ; incorporating state-of-the-an 
computational fcatSes of tire ^ codTind^“ “L S g U ! ‘ apen!d and “W* «***• R<*>vam 
general Body Fitted CoordSSewrtS “£L 21 ^ y/n ? tatmg coordinates, cylindrical and 

variable atnigement ™dvS ESTrlu 8 - 0Kler “hemes, colocated 

moving grids advanced turbulence models, incompressible/compressible flows, and 

iSSc^reS development, code demonstration for 

for genetic aniSSS.^S\7S^K n m S?*- *£' “““ loss coefficients 
configurations. P extensions to labyrinth, damping, and other seal 


INTRODUCTION 

tot^gKStlS re fiS SJjSfia* ^Progressively higher efficiencies, 
the most compleT enrine coraSnent JS, fi ? ^ a lo ? SCT Ufe • Turbomachinery, 

leakage, rotodynamic respome Sch^c^cimnnnlnf^^ Seal con ^g ur ?tions to control 
incompatible fluids. As all of the ahnve aff»rve p - en( lQ icrances, and to inhibit mixing of 

or limit cycles, it is important to understand the^d^mS^Tfl 6, thro ^ gh P arasitic tosses 
their interdependence As modern flow 311(1 rotodynamics and 

increasingly fiX? S^thffaS^ , becomes lighter and is operated at 

important deS|n Jon?em Prime , SS£/^ t t u aa • thc rotod y namic stability becomes an 
Shutfle Main Engine (SSME? Fuel S ° f , seals are in th e Space 
becomes a critical element for the staS nation n tarbo P um P s where seals perfonnance 
efforts of NASA LeRC, MSFC and Texa« ajJm rr f engin ® s * Coordinated seals research 
Shuttle's reliability £d perfo, ro^ce hfs haVC to the 

and seals rotodynamic coefficients reauires advan^.^ 1 accurate . prediction of leakage rates 
complexity and extreme computational models. Geometrical 

and bulk flow models to theirSts In 1988 tJ^SISS 011 fll w s pushed existi ng empirical 
and Mullen 1 demonstrated Sat CRD cm Muszynska - Hendricks, Brown 

and can complement the analysis of rotodynamic instability oK£ mpIeX fl ° W 111,1011811 
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The basic objectives of this program follow. 

1 To develop a verified three-dimensional “scientific” CFD code for analyzing 
seals, representing the state-of-the-art in accuracy of predicting flow fields 
rotodynamics and performance of a wide variety of seal configurations. The 
code will provide insight into flow and dynamics, contribute to the technical 
data base and provide an accurate standard for less complex industrial 

codes * 

2 To compile and generate a set of verified one-, two-, and simplified three- 

dimensional “industrial” codes that will enable the designer or researcher to 
expeditiously analyze the wide variety of seals. 

3. To develop a modem seal design computational environment with 8 

Based System (KBS) and links to CAD/CAM packages and scientific data 

visualization tools. 

The protect team members include NASA LeRC, Mechanical Technology Incorporated 
developing industrial codes and KBS, and CFD Research Corporation developing a scientific 
CFD seals 8 code A parallel research effort is being conducted in which the scientific code will 
contribute to the technical data base and provide an accuracy standard for less complex 
codes. The industrial codes will provide speedy analysis of seal designs. Graphical User 
Interface (GUI) and KBS data base will be a unifying mechanism. The technology tr ^s c 
will be accomplished via annual workshops at which the program workplan will be reviewed 

and revised. 

The aim of the code development effort is to develop a steady/transient 3D CFD code for the 
analysis of fluid flow for a broad range of seal configurations and for the ^ ial >' sl . s ^ t J 1 . e 
dynamic forces on the seals. The geometrical seal configuration will include cylindrical, 
labyrinth, honeycomb, damper, brush, face, wave, grooved, tip, contact, and noncontiguous 
seals The fluids will range from incompressible liquid, two-phase mixtures, to compressible 
gases An existing 3D CFD code, REFLEQS 2 , is being adapted for seals by incorporating 
physical models to test turbulence, non-isotropic wall roughness, phase change, and 
Sation The scientific code was demonstrated at the 1992 annual workshop m August 

1992 at NASA LeRC. 

CFD CODE ENVIRONMENT 

Existing CFD methodology of turbomachinery seals design is an expensive and time 
consuming process. Geometry definition, grid generation, selection of boundary conditions, 
parameters selector., post-processing, and analysis of results often require CFD 
expert involvement. Rapidly advancing computer system software s ' a i 

allow the automation of the CFD-based design process. In order to develop CFD-based seal 
ddp“ ent integration of CAD/CAM, CFD, structural mechanics, rotodynamics. 
computer graphics under lie Knowledge Based System (KBS) is currently conducted. 

The tasks of KBS are to integrate scientific and industrial codes with a consistent database, to 
provide a user-friendly graphical interface with context-sensitive help and to provide the 
Expert System to guide the seal design and optimization process. A major component of the 
KBS environment 8 will be the scientific CFD code. The outline of the structure and brief 
description of the components of CFD seals code are provided. Technical specifications of 
the CFD code capabilities are discussed in the next section. 

Figure 1 presents a flow chart for the CFD code design environment 

svstem is the CFD code and most of the current development effort is dedicated to it. The 
oUicr major modules CFD GUI -Preprocessor and KBS/GUI are being developed in 
collaboration with Mechanical Technology, Inc. The remaining components, CAD and 
Visualization will be available from commercial vendors and the convenient links to them will 
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provided within GU?-Preprowssor°^ 30(1 aIge .^ ra ' c gr id generation will be 

such as ICEM-CFD, PATR AN, or IDEAS will C f Dpackagc 

(SdsTiii'mmD Jnnf?/ Cslablishcd intcrf accs with PATR AN and ICEM S -CFD 8 codcs l 

standards will be adopted *i^rfaSl f‘ rcc f ,y " ilho . ul any translation. Similar 
format files for grid vectors J,d S, scientific visualization tools. PLOT3D 

sipercomputeiT “ y UNIX X - w,ndow based workstation or on mainframe or mini- 
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Figure 1. 


KBS - ^.edge Base 

CFD CODE CAPABILITIES 

™ .^SdeSlS’flow Womatonfora ^, in ! ius “ ial ***■ codes, 

die details of flow tmd i rerfSES J VJ -fF ° f ?*■ “ w,u 3)30 be used to predict 
inadequate testing has been performed With . con ^Sy ratlons as well as seals where 
following features as a minimum: d ‘ W ^ ^ m m,nd ’ a sc,entific code must possess 

StSi^ a a V rJ r 1nH 0ke a S in generalized body-fitted coordinate grids- 

Stationary and rotating coordinate frames; 8 ' 

Steady-state and time-accurate solution capability- 

incompressible and compressible flow solutions- 

cSo„ P S; ProPertieS (ViSC ° Si,y - *“'*• ^ >**< «C )i 

Provision for stepped surfaces and injection ports; 

Inclusion of viscus dissipation in energy equation; 

Treatment of sources due to external fields; 

Variable surface roughness treatment; 

Provision for effects of pre-swirl and upstream effects' 

Customized input and output features for seals; and 
calculation of rotordynamic coefficients; 

was selected. ^This codeT^^ *J° W E 9 uati on Solver) code 

grants as well as IR&D funds. Severn! 


1. 

2 . 

3. 

4. 

5. 

6 . 
7. 

9. 

10 . 
11 . 
12 . 

13 . 

14. 
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Some of the features of the code arc listed below. 

1 High-order spatial differencing schemes including central-differei 

several third order schemes such as Osher-Chakravarthy, Osher, Superbee, 

2. Sec oiid -order temporal differencing: Crank-Nicholson and three-point 

3 rS~s-f boundary conditions: specified mass flux, specified 

pressures waU and symmetry. Some seal related conditions include specified 
K pressure condition wiih inlet loss factor, rotating and whirling wall 

■ 32Effi82S5**ft=“ - " 

assumed to undergo whirl in a circu a mi a<d stead v The fluid pressures on the 

frame whirling with the rotor renders me P^.-^^^orine forces at several different whirl 

SsiSSHB-JKSs — 

eccentricity as well as rotor misalignment. 

SELECTED RESULTS 

validation of the code, and to ensure code reliability. These are. 


1 . checkout problems; 

2. benchmark problems; 

3. validation problems; and 

4. field problems. 
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Sd fnsurelhe' accu^v *Sr ,tf"„ , pul th ™ gh J ln “*«*« army of lest problems to assess 

Fully-developed and developing now in pipes, channels, and long narrow 

laminar now between rotating cylinders, with and without Taylor vorticess 
r™,e„. T cav “)' n”'*': (Re up to 10,0007); 3D driven cavity now ' 

bearing*” WlUl 2nd Wlth0 “‘ heat ,ransfer Planar wed|e How in a slider 

Laminar flow over a backstep*. 

Laminar flow in a square duct with a 90° bend? 

Turbulent flow in a plane channel i°. 

Turbulent flow induced by rotating disk in a cavity* 11 
Centripetal flow in a stator-rotor configuration* >2 
Turbulent flow in an annular seal 1 *. 

Turbulent flow in a 7-cavity labyrinth seal 14 . 

Tutbulent compressible flow and heat transfer in turbine disk cavities 1 * 

Wilh lid clearance “ d «W prossuro Siera and 
control of flow through vortex imposition 1 **. 8 ent ana 

p mfinite ^finite length bearings (without cavitation) 1 ?. 1 * 

?omn d sS C ,. C ^i!i C l e ?,'J a lc “ lati ° ns for incompressible flow sealsnuo and 

Computations performed under NASA MSFC Contract No. NAS8-38101) 


2 . 

3. 

4. 

5. 

6 . 

7. 

8 . 

9. 

10 . 
11 . 
12 . 

13. 

14. 

15. 

( * 


P* = Pi 


in 


- 1/2 (1 + £ ) pu.2 


order to estimate these ncinn »h 0 ^ ^ senes of cases were computed in 

conditions. These results are also defcribcd P £tow ‘ ,iffl:rcnt n ° w and JCbnrotry 

Swero’SrodC M.5 ™ ‘Sw SS m'“ n , The Xal P a ™"®* 

mm, inlet pressure = 1.52 MPa sump orcssure M n ' n ^i an<l ex,t = 0172 and 0 086 

=1.4, dynamic viscositv — i s y inT pff! L e ~ 065 Mpa ’ mlet temperature = 650 K, v 
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and the computational grid used had 12, 6, and 30 cells in the axial radial and 
circumferential directions. The results for this case are shown m Table 1. Also shown arc the 
results from bulk flow theory as reported in Reference 21. As seen, the two sets compare 
reasonably well, including the prediction of flow choking at the exit of the short seal. 


Tabic 1. Tapered Gas Seal Rotordynainic Coefficients 


UD 

KN/m 

k N/m 

CN-s/m 

c N-s/m 

Exit Mach 
Number 

0.1 

948000 

19700 

11.5 

0.016 

\ 1.00 ' : 

1150000 

15429 

9.94 

0.043 

1.00 

A 9 

1700000 

75700 

43.9 

0.073 

0.96 


2125500 

60886 

38.62 

0.09 

0.97 


2880000 

267000 

152 

0.27 

0.83 

0.4 

3553200 
— tt— yr 

233820 

145.9 

0.57 

Iac*i iHr 

0.83 

L "" PA-92-24H 


1 Nelson, 1985^ | I Present Results 


Incompressible long seal: The seal geometry for this case, as described in Reference 20 is : 
seal radius = 39.656 mm, clearance = 0.394 mm, length = 240 mm^ Properties of ^ 
20°C were assumed: dynamic viscosity = 1.04 x 10' 3 Pa-s, density - 998 Kg/m . g 
used had 20, 15, and 30 cells in the axial, radial, and circumferential directions, respectively. 
The low Reynolds number k-e turbulence model was used in the calculations. Thewhiriing 
rotor method was used to calculate the rotordynamic coefficients. The calculations were 
performed at three shaft speeds: 1080, 1980 and 3000 rpm. Six pressure differenuals fo 
each speed were used: 20. 50, 10, 250, 500 and 900 KPa. Inlet swui was specified 1 using the 
experimental data. Results of these calculations are plotted in Figures 2 through 6. 
Experimental data from Reference 20 are also plotted on these figures for comparison. TTie 
comparison between the experimental and computed values is very good; the negative direct 
stiffness coefficients at high rpm and/or low pressure differentials are also predicted correctly. 
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Figure 4. Direct Damping Coefficient. Figure 5. Cross-Coupled Damping 

Symbols represent experimental Coefficient. Symbols represent 

data form Reference 20. experimental data from 

Reference 20. 



Figure 6. Direct Mass Coefficient. Symbols represent experimental data from Reference 20. 


Entrance loss coefficient: As described above, the entrance loss coefficients is an empirical 
parameter which depends on several factors such as seal entrance geometry, flow rate etc To 
analyze the effects of some of these parameters on the loss coefficient, a limited parametric 
study was performed. A generic annular seal geometry was assumed, as shown in Figure 7. 
The flow was assumed axisymmctric. Effects of three parameters was considered: 1) radius- 
to-clearance ratio, 2) entrance region radial gap to clearance ratio, and 3) flow rate. The 
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nominal values were: seal radius = 25 mm, seal length = 25 mm, and the axial length of the 
inlet region kept at 1 .5 times the seal radius. 


Entrance Seal 



Figure 7. Flow Domain and One of the Grids Used for the Entrance Loss Coefficient 
Calculations 


Water at 20°C was used as the working fluid, and the rotational speed was kept at 3000 rpm, 
•which corresponds to a rotational Reynolds number of 1.88 x 10®. Three radius to clearance 
ratios were considered: 50, 100, and 150. Two gap-to-clearance ratios were used: 50 and 
100. Each of these combinations were run at several different flow rates (axial reynolds 
numbers). The axial and rotational Reynolds numbers were defined as: 


and 


Re a x — 2 ■ (CL)u w>m /v 
Reg =G)r 2 /v 


where v is the kinematic viscosity, u M m is the mean axial velocity, CL is the seal clearance, r 
is the seal radius, and to is the shaft speed in rad/s. The exit pressure of the seal was kept Fixed 
at 0 Pa, and nominal flow rates were imposed at the inlet. A fully-developed turbulent axial 
velocity profile was imposed at the seal inlet. Computational grids used had 5 cells in the seal 
clearance, and 30 or 50 cells in the radial directions in the entrance region depending on the 
radial gap to clearance ratio. In the axial direction 30 cells were used in the entrance region 
and 20 cells in the seal region. Appropriate clustering was used to resolve the critical flow 
regions. 

To calculate the loss coefficient, the seal entrance pressure is required. Moreover, this 
pressure should correspond to the value which will generate the specified flow rate in the seal 
when only the seal is considered. To estimate this, pressures in the downstream part of the 
seal length were used. The flow becomes fully developed in this region, and the p ressu r e 
varies linearly with axial length. This slope was calculated and extrapolated to the seal inlet to 
provide the inlet static pressure. The loss coefficient was then calculated based on this 
pressure, and the computed pressure at the entrance of the flow domain. 


r 
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results at compile axT^I^ ggf ° f °* 

2. ,ncrcasing 

radius ratio, but the sensitivity is lower g er entrance region gap to clearance to 

3. for a fixed geomeby, the loss coefficient reduces with higlKr „ ow ralcs 

For die configurations tested, die loss efficients were between 0.406 to 068. 

.Table 2. Entrance Loss Coefficients, Ra dius/rien„.w = 30 

Ifift fnanyTln-5 PA T 11 — 

Entrance Gap/Cleara nce = 100 

I r~ — i — 


u ax 

R %x 

% 

10.814 

16.232 

21.619 

26.942 

10377 
15484 
1 20746 

25854 

0.471 

0.431 

0.414 

0.406 


u ax nVs 

Re ax 

% 

10.82 

10.24 

21.66 

27.06 

10384 

15584 

20785 

25970 

0.490 
I 0.488 

0.482 
0.48 


PA-02-24 Q 


. ” ^micuiLC JUUSS coerii 

1 Entrance Gao/CiMr-anAA - ca 

cients, Radius/Clearance = 100 

"Vs 

Re ax 


I — Entrance GaD/ClearancA - inn I 

% 

Uax mfs 

Reax 

5 

10.80 

16.56 

21.595 

26.67 

32.27 

43.062 

5181 

7945 

10361 

12796 

15484 

20667 

0.562 

0.54 

0.526 

0.51 

0,493 

0.478 

10.797 

16.176 

21.55 

26.934 

32.24 

42.533 

5167 

7761 

10339 

12664 

15469 

20408 

0.567 

0.558 

0.55 

0.54 

0.537 

0.524 


Entram 

Uax 01/8 

i A QA 

:e Gap/Clearanc 
R %x 

VrVArUi 

e = 50 
\ 

Kaaiuj 

Entram 

Uax 

s/Clearance = 15 
ce Gap/Clearanc 

Re ax 

0 

e- 100 
% 

* U.Oc 

16.19 

21.49 

26.74 

32.25 

48.33 

64.487 

3461 
5178 
6874 
' 8553 

10315 
15461 
20630 

0.66 

0.65 

0.647 

0.637 

0.628 

0.606 

0.595 

10.75 

16.09 

21.47 

26.81 

32.176 

47.87 

64.165 

3438 
5146 
6874 
I 8553 

10292 
15315 
20630 

0.68 
0.66 
0.65 
! 0.648 

0.64 
0.63 
0.624 


PA-8Z-24 M 
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CONCLUSION 

Advanced, CFD based, tu rboma ^‘" Yncludfne ^sdenUfic CTO^ow solver, and a range of 
consists of several major compoi *™ lud A ical user interface will be integrated with a 
industrial seal design codes u t0 C AD/CAM and 

Knowledge : B ased System i and mU of inte g rated CFD des *f 

scientific visualizabon tools. The pape p ™ w solver as major component. The 

code system under the KBS with m etho& has been 

scientific CFD flow solver incorporating state-of-^art ™™ e £ c monstrated for seV eral 
developed, validated on a broad range , p hag also been use d to calculate the 

cylindrical seal flow configurations. ocmoared with experimental and other 

rotordynamic coefficient of number of t study^show very good agreement with 

analytical/numerical results. Results of valida ^/^^J^eLremfnts. A short 

ssss.ss&’S.' sas. .0 

J-SSra epee gap to sea. Ceatance and me Hoe, ran, 

also change the loss coefficients, but to a smaller extent. 

The present efforts include cxten d ing t^e CF^ codc^^o ilab^nnth^d^dam^r^sca [Q 

c^er^^broad^ range °of* seal S «mfigutati 0 ns for fluids ranging from incompressible to 
compressibel gases in two-phase flow. 
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